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Abstract 

Under the steady state condition, the spectrum of electrons is investigated 
by solving the continuity equation under the complex radiation of both the syn¬ 
chrotron and Compton processes. The resulted GRB spectrum is a broken power 
law in both the fast and slow cooling phases. On the basis of this electron 
spectrum, the spectral indices of the Band function in four different phases are 
presented. In the complex radiation frame, the detail investigation on physical 
parameters reveals that both the reverse shock photosphere model and the for¬ 
ward shock with strong coupling model can answer the a ~ — 1 problem. A 
possible marginal to fast cooling phase transition in GRB 080916C is discussed. 

The time resolved spectra in different pulses of GRB 100724B, GRB 100826A and 
GRB 130606B are investigated. We found that the flux is proportional to the 
peak energy in almost all pulses. The phases for different pulses are determined 
according to the spectral index revolution. We found the strong correlations 
between spectral indices and the peak energy in GRB 100826A, which can be ex¬ 
plained by the Compton effect in the fast cooling phase. However, the complex 
scenario predicts a steeper index for the injected electrons, which challenges the 
acceleration mechanism in GRBs. 
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Subject headings: gamma-ray burst: Band-like spectra: Klein-Nishina regime: 
Synchrotron radiation and Compton scattering 


Introduction 


Althongh a remarkable advance of investigations on GRBs was made in the past decades, 
the nature of the prompt emission of GRB is still unclear. Most spectra of GRBs can be de- 
scrib ed by the empirical Band fnnction, whose format is a smoothl y connected broken power 


law fiBand et ahl Il993l: iGoldstein et al.l l2012l: iGruber et ahl 120141) . Many ph ysical models 


have been built to explain the radiation mechanism in the prompt phase. iBeloborodov 


(120lOl) discussed the thermal photosphere emission from the passively cooling jets. Inelastic 
nnclear collisions pl ay an important role in producing the mnltiplicity of high energy photons 
(iBeloborodovI 120031) . The resulted spectrum roughly is a broken power l aw except a high 


energy cut off, and the range of the low energy spectral index is limited. iRees fc Meszaros 


(120051) considered the leptonic dissipation model to acconnt for the high Inminosity of GRBs. 
One advantage of the photosphere model is the prediction of the observed peak photon en¬ 
ergies. However, the detailed spectral indices can not be read directly from the thermal 
radiation to fit observations. By including the reasonable ingredients like the bnlk motion 
of the jet, the location of the dissipation and the viewin g angles, the phot osphere emission 
can also lead to a variety of non-thermal spectral shapes (IRvde et al.ll2nil). The more pop- 
nlar i dea abont prornpt ern ission of GRB is the relativistic fireball model (IRees fc Meszaros 
19921: iMeszarosI 120021120061 ). An extreme fast ontflow collides with the interstellar medinm 
and produces a violent shock. Electrons in the plasma will be accelerated by the shock 
to form a power-law energy distribntion via Fermi acceleration (jPeacockI Il98ll ). The syn¬ 
chrotron radiation (SR) of these electrons will lead to a broken power law spectrnm in the 
op tically thin reg i me. T he magnetic reconnection and turbulence model has been proposed 
by IZhang fc YanI (120111 ) to explain the GRB emission. This model can overcome the low 
efficiency and electron excess problems of the internal shock model. Thanks to the great 
progress on the detection of high energy astrophysical photons, especially the Fermi and 
Swift missions, these models can be well testified by the spectral analyses of GRBs. 

The Fermi Gamma-Ray Burst Monitor (GBM) spectral catalog of the first two years 
and fo ur years indicated that the broken power law is still good in fitting most spectra of 
GRBs (IGoldstein et al.ll2012l: iGruber et al.l 120141) . Here, we denote the low and high energy 
spectral indices as a and /3, respectively. For the high quality “BEST” samples, the low 
energy spectral indices a peak at —1; Up to 17% samples violates the synchrotr on —2/3 “line 
of death”; while additional 18% exceed the —3/2 synchrotron cooling limit (IGruber et ah 
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2014J) . This challenges the standard shock plus sync hrotron models. Other factors need 
to be considered to overcome or recoil this problem. iTvnndman et al.l (120131) argued that 
jets with angle dependent Lorentz factor prohles can produce a ~ — 1, where photosphere 
emission of non-dissipative jets are considered. This suggests us that the thermal emission 
plus the geometry effect can create the observed GRB spectra. Another remarkable model 
which can solve t he a —1 problem is the fast cooling synchrotron radiation in the decaying 


magnetic field (jUhm fc Zhandl2014l) . The low energy spectral slope a is found t o be related 


to the decaying slope of the magnetic filed by an analytic method. Further on, I Zhao et ah 


fl2014j) presented the GRB spectrum with the decaying magnetic filed by considered both 
the synchrotron and inverse Gompton (IG) radiation. The IG component is also important 
in generating the low energy spectral index. iDuran et al.l (120121) argued that the inclusion of 
the IG component will harden the spectrum. But, it is still difficult for a to approaches the 
observed value —1, and the maximal value of a is obtained in an extreme condition. In this 
paper, we will show that the a ~ — 1 problem can be solved in the complex of synchrotron 
and Gompton processes, where —1 is not a limit but an median value of the possible range. 
We also investigate the time resolved spectra of several bursts, and show the evidence of the 
complex radiation mechanism. 


Nakar et al.l (120091) discussed the synchrotron and synchrotron self-Gompton (SSG) spec¬ 


tra with the Klein-Nishina (KN) effects, they showed that the spectral energy distribution of 
electrons becomes hard, and the observed synchrotron spectrum is hardened. Our work here 
solves the continuity equation analytically, and the spectral indices of electrons are obtained 
directly. Then, we present the full GRB spectra in the complex situation. We show that the 
solution of the a ~ — 1 problem admits quite physical parameters. The non-trivial relations 
between indices and peak energies will be presented, which can be considered as a sign of the 
complex radiation. The organization of the paper is arranged as follows. In Section 2, we 
investigate the electron s pectral distribution in the complex of synchrotron and IG, which 
completes the analysis of iDuran et al.l (120121) . In Section 3, we derive the corresponding 
GRB spectra in both the fast and slow cooling cases. In Section 4, we present all possible 
radiation phases and corresponding spectral indices. The analysis of time resolved spectra 
is given in Section 5. Discussions and conclusions are given in Section 6. 


2. Continuity equation of electrons 

Most works investigate the photon spectra in the prompt phase with the assumption 
that the injected electrons have a power-law distribution and cool down instantaneously. The 
observed emissions probably originate from electrons in different emission states and regions 
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flDuran et al.ll2012h . After interacting or radiating photons, the spectrum of electrons in the 
plasma changes. The electron evolution is governed by the continuity equation, which is an 
reduced version of the Fokker-Planck equation for electrons. The ignored terms are related 
to the dispersion and escaping effects. The Larmor radius of electrons is much smaller than 
the emission re gion, and the disper sion effect is not important in the GRB scenario. In this 
work, we follow iDuran et al.l (1201211 to solve the continuity equation analytically. In the local 
jet fram^, the continuity equation is written as 


4]v(y) + A [yA,(.,')] = s(Y), 


( 1 ) 


where 7 ' denotes the Lorentz factor of electrons. 5 '( 7 ') denotes the energy distribution of 
the source electrons. It is widely accepted that the shock-accelerated electrons have a power 
law distribution, i.e.. 


S{^') oc 


7 


/-p 


v' < v' < 'y' 

/m — / — /m 


7 ' < 7m, 


or 


7'>7m 


( 2 ) 


For shock accelerated electrons, 7 ^ = ^T (F is the bulk Lorentz factor of the je t, g, 


is the constant fraction of the shock energy) is the minimal Lorentz factor of electrons fjSari 
199811 . The maximum energy for electrons gained in the sh ock acceleration via Fermi process 


7 h.,.^ is written as 7 ^^,^ ~ 4 x lO’^il' f Kumarl let ah 20121: Fan k. Piran 


200811 . 


For a steady case, one sets N{'y') = 0. The electron will lose energy via synchrotron 
radiation and IC scattering processes. The cooling of the electron is described by 


- 7 ' = 




+ 


O'KNL'f 


n 




( 3 ) 


where (Tkn is the cross section valid in the KN range. Dehning rj = y'/w^seed/dTeC^ (t'seed is 
the frequency of the seed photon), rj describes how deep an electron is in the K N regime. In 
the T homson regime, one has 7 <A 1, while one has 77 > 1 in the K N regime flFan &: Piran 
2008 1. Setting axN = 07 /( 7 ), /(h) is written as ( Duran et ah 2012 1 


/(>,) = j 


l + r] f 27(1 7 ) 

y \ 1 -|- 27 


ln(l + 27])] + 


ln(l -|- 27]) I + 3r] 


27] 


(1 + 27)2 


( 4 ) 


However, one electron loses energy ~ in the KN regime, and the power of 

Compton scattering is proportional to 7 ', rather than y"^ in the second term of Equation ([3]) 


^Here we use the prime sign to mark parameters in the jet frame. 
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flFan et alJl2005h. We need to presen t a formula which is valid both in the Thomson and KN 
limit. Following iFan fc PiranI (120081) . we dehne a f{r]) function, i.e., 


fiv) 


fiv) 

l + r]' 


The general Compton-Y parameter is dehned as 




fiv) 


U'b' 


(5) 

( 6 ) 


where U'^ = L ^/and U'^ = refer to the photon and the magnetic held energy 

density, respectively. Yq describes the ratio between the power of the Compton process P{q 
and that of the synchrotron radiation in the local frame. Now the cooling rate of the 
electron is rewritten as 

- 7' = ^(1 + Yc). (7) 

-^syn 

where T'y^ = GirMec/aTB'"^. In the complex of the synchrotron and IC processes, electrons 
will lose energy signihcantly if the cooling time is smaller than the dynamical time ^ 
P/Tc. When they are the same, the critical Lorentz factor of electrons is expressed as 


7c 


Q-kYMpC^ 


arPB'^ 1 + Yc 


( 8 ) 


c 


Duran et al.l (120121 1 mainly investigated the spectral slope in the low energy range. Thus, 
the electron spectrum was solved analytically by the continuity equation in the low energy 
range. We will complete the analytical solutions of the continuity equation in the full energy 
range. We also consider bot h the fast and slow cooling pha ses, analogous to the full spectrum 
of the synchrotron radition (jSarilll998l: iFan fc PiranI 120081) . 


Considering the stationary case, one sets ^(y') = 0. The continuity equation is written 


as 


d'y 


7 [7'iV(7')] = S(7'). 


(9) 


The analytical solution is given by iV( 7 ') cx 7 ' 


7 ' < 7 ' < 7 ^, one has N 
energy range is given by (IDuran et al.l 120121) 


J S{'y')d'y'. In the fast cooling phase, 
7 O PC 7 ^ The c orresponding spectral index of electrons in this 


d IniV(y) ^ dlnfjr]) Yq 
d Iny' d Iny 1 + Yc 


( 10 ) 
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The formula of d\af {rj)/d\n.r] can be approximated in the two limits, i.e., 

d\rif{ri) _ r -3?? h < 1, 

d\n.r] 1—2 S> 1 . 

One observes that the value of (iln/(? 7 )/(iln ?7 is always negative. Since the Compton param¬ 
eter Yc is always positive, the range of pi is 0 < pi < 2. This result have a signihcant impact 
on the spectral indices of GRBs. 


In the slow cooling phase, i.e., 7 ^ < 7 ' < 7 ', the electrons lose no significant energy 
via radiation. So, the spectral distribution of electrons is roughly invariant. This leads 
to N{'y') oc 7 '“^. The distribution of electrons in the high energy regime is also changed 
according to the continuity equation. For 7 ' > max( 7 ^, 7 '), since the source has the form 
S'( 7 ') oc 7 '“^, the corresponding power-law index can be obtained analytically 


dlniV(y) ^ ^ dlnf{7]) Yg 

d Iny' ^ d Int] 1 -|- Tc 


( 12 ) 


Following the argument previously, the ra nge of n? is (p — 1) < p 2 < (p -|- 1). Therefore, our 
work here completes the analysis given by iDuran et ahl (120121 ). Collecting these results, the 
distribution of electrons is written as 


A^(7') oc 


j'-Pi 7 ' < 7 ' < 7 (^, 

7'“^ 7m < 7 ' < 7c, 

max{o;„, 73 < y < 7i,„. 


(13) 


The cutoff factor 7 ^^^^^^ is not analyzed here, since it occupies a negligible fraction. Setting 
Yc = 0, which means that the Compton processes are ignored. Equation (1T^ reduces to the 
synchrotron radiation case. From Equation ([13]), it is evident that the inclusion of Compton 
processes changes the spectra of electrons, and further affect the observed photon spectra. 


3. The spectra of photons 


The photon spectra in the prompt phase of GRB can be studied on the base of the 
electron spectrum in Equation flT^ . First, we discuss the low energy synchrotron spectrum. 
In the fast cooling phase, the distribution of electrons is composed of two connected power 
laws with indices —pi and —p 2 , respectively. It is well known that the flux of synchrotron 
radiation for the low energy range is proportional to One way to estirn ate the flux in 
the low energy regime due to the synchrotron is given bv iDnran et al.l (120121 ). i.e., 

,, H 1/3 


POO 

Fy“ = ^(R') / d7'iV(7') 

J'y', 


Z/(7') 


( 14 ) 
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where 7 ^ = ^/27rM^a7Je&, and jy{'y') is the frequency corresponding to 7 '. The integration 
fT is done for 7 ' > 7 ^. If 7 ^ is smaller than 7 ' then the integration is composed of tree 


parts, i.e., Jy + The result of the integration is independent of 7 ^, since there 

is no electrons distributing for 7 ' < 7 '. Thus, one always has Th oc IPnran et ahl ( 20121 ) 
considered that there is no lower cutoff for the electron spectrum, and obtain a formula of 
the spectral index. In the presence of bounded electron spectrum, their formula is valid only 
for 7 '. We present the formula of the spectral index at this critical frequency, i.e.. 


dlnF,Y° ^ 1 _ (5/3 + q)/2 _ 

dlnuc 3 l-(p-l)/(p + 2/3 + g)(7/y7(„)5/3+T 


(15) 


where we define q = for concision. The formula here agrees with Duran et ah 


dlnrj l+Vc 

(120121 ) except that we take use of the formula of p 2 in Equation (IT^ . etc is meaningful to 
describing the spectral curvature. But etc can not describe the full spectral slope at the 
turning point, since only the low energy tail of the synchrotron emission is considered to 
obtain et. 


When 7 ' < 7 / < 7 (,^, the integration contributes to the spectral index. However, 
Equation (ITT)) counts only a small part of the synchrotron radiation power. The photon 
spect rum should be derived via the relation F^du oc N{'y')P'dj' and z/ oc 7 '^ (jFan fc Piran 


20081 1. By considering these relations and combining the synchrotron self-absorption, the 


full photon spectrum in the fast cooling phase is written as 


Ff ° oc 


V 

1 

1/3 

PI -1 

V 2 
P 2-1 

V 2 


y < 

Z/a < Z/ < Z/c, 

Z/c < Z/ < Z/m, 

^ ^ ^max* 


( 16 ) 


Ua refers to the absorption cut off frequency. The presence of Compton process changes 
the second index of the spectrum, and this breaks the “line of death” of the index of the 
synchrotron radiation. 

Now we discuss the spectrum in the slow Icooling case. According to Equation (fT3|) . the 
spectral indices of electrons are —p and —p 2 for the low and high energy range, respectively. 
Therefore, the synchrotron induced spectrum of photons is written as 


F/y" oc 


Z/3 

P-1 

z/ 2 

P 2-1 
U 2 


y < t'a, 

< Z/m, 

Prn<y < yc, 
yc<y < i^max- 


(17) 

















Compared to Equation ffT^ . we only replace pi with p in the slow cooling phase. The typical 
value of p is larger than 2 for the shock accelerated electrons. We also have 0 < pi < 2, see 
last section. Thus, the fast cooling case has a harder spectral index than the slow cooling case 
in the intermediate energy range. For the highest energy range, the spectrum becomes soft. 
The reason is that the most energetic electrons mainly lose energy via Compton processes, 
and their synchrotron emission is not as efficient as the pure synchrotron radiation case. 


The Compton process will produce a spectrum in the energy range above the synchrotron 
spectrum. The source photons of the Compton processes can be of the synchrotron origin. 
In this case, our scenario mainly discusses the synchrotron self-Compton (SSC) process. We 
estimate the energy of scattered photons in order to understand whether these photons can 
be observed in GRBs. The typical bulk Lorentz factor of GRB jets in the prompt phase is 
several hundreds. The observed peak energy of most GRBs is about hundred keV. In the 
local jet frame, the energy of synchrotron photons peaks at the order of keV. Considering 
the moderate shock acceleration efficiency e ~ 0.1, and p = 2.5, one obtains ~ 6000. 
The T] parameter is calculated to be r; ~ 12. So the SSC process happens in the KN regime. 
The scattering process will transit half of the electron energy to the photon. Thus, the 
observed photon energy can be as high as 150 GeV, which can be recorded by the LAT 
monitor. However, such high energy photons can be attenuated by three significant effects. 
The first one is due to the cross section in the KN range, which strongly suppresses the 
collision chance. S econdly, the optical depth for such high energy photons can be l arge in 


the prompt phase f Bosniak fc Kumar 2012 : Beloborodov 2010; Chang et ah 2012^. T his 


will lead to the time lag phe nomenons for GeV photons in GRB Abdo et ah f 2009allbl) : 


Ackermann et al.l fl2010l. 120111 1. The production of such photon collision is pairs of electron 


and position. This may enlarge the radiation populations in the jet slightly. The third effect 
is caused by the extragalactic background lights (EBL). The absorp tion can be described by 
the model dependent gamma-ray opacity flArchambault et al.ll2014l l. This effect suggests us 
that GRBs can enlarge the positron population in the universe. All these effects reduce the 
chance of VHE photons to arrive the detector. 


The power of the Gompton processes in the GRB is about the same as the synchrotron, as 
indicated in the next section. Although the cross section suppresses the collision probability, 
the Gompton scattering can not be ignored only by this reason. The Gompton induced 
spectrum is most evident in blazars, which is also produced the radiation i n jets of AGNs . 
The spectral energy distribution (SED) of Blazars shows a bimodal pattern (IBonnolillet ahl) . 
The first peak in SED, usually in the optical to the soft X-ray range, is due to the synchrotron 
radiation, while the second peak in the high energy range is of the SSG or external Gompton 
radiation. If electrons in the jet have a complex spectrum as in Equation flT^ . one can 
expect that the spectral indices of photons are strongly related. Since we concern the GRB 
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spectr um in this work, the s tudy of the complex spectrum for blazars will be given in another 


work fl Jiang et ah. Ilin prepf) 


4. GRB spectra 


Most GRB spectra are described successfully by the empirical Band function. Now, we 
aim to obtain the spectral indices of the Band function by using the Fy spectra in Equations 
ffTbjl and ffT7|) . The empirical Band function is a smoothly jointed broken power law, which 
reads 


(a+‘2)E 


N{E) = Ax 


lOOkeV 
■ E 
lOOkeV 


)/3e/3-a 


( j 

VlOOkeV ) 


OL — ^ 


E < Ec, 

E > Er.. 


(18) 


keV ^ flYu et ahll2014l) . 


A is the normalization factor at 100 keV in units of photons s“^ cm“^ 

The relation of Ep and is Ep = Ep{a — (3)/{2 + a). Some bursts show a double broken 
power law spectra. In orde r to analyze su ch anomalous spectra, we introduce the third 
spectral index 7 . Following mi et al.l (120141) . the empirical function for the double broken 
power law is written as 


1 

f (f )“> 

E < Ehi, 


N{E) = Ax{ 

1 ( Ebi\a-0( E '1/3 

1 E„ ■> '' En ^ ’ 

Ebi < E < Eb2- 

(19) 

1 


E > Eb2- 



Here A is the normalization flux at certain energy 


By the relation E'^N{E) one obtains the expression of a, (3 and 7 as functions 

of Pi and P 2 in both the fast and the slow cooling phase. There is still another phase that 
Vc and Vm is too close to distinguish, which is named as the “marginal” phase or the “both” 
phase (IYu et al.ll2014j) . Since pi and p 2 are functions of p and g, one can express p and q in 
terms of a and (3 for a given phase. In this way, p and q can be obtained from observations. 


For the Band function, one remaining important question is that Ep corresponds to 
which frequency, Um or z/c. Since Ep is the peak energy of the uF^, spectrum, one concrete 
criteria is that the spectral index of uF^, spectrum is positive for E < Ep and negative for 
E > Ep. Note that this criteria is not valid for the double broken power law. In the fast 
cooling phase, Ep is equal to hvrn) because 0 < pi < 2. Meanwhile, p 2 should be larger than 
3. This can be used as a consistent condition to constrain the physical parameters. In the 
slow cooling phase, Ep is produced by electrons with 7 („ for p > 3. If p < 3, Fp corresponds 
to hvc- So, there are two possible scenarios to explain the Band function in the slow cooling 
phase. A double broken power law is possible in both the fast and slow cooling phases. Both 
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Ebi and Eb 2 are uniquely determined if the absorption part is not included in the spectra. 
We will not discuss the synchrotron absorption effect unless some extreme hard low energy 
spectrum needs explanation. Collecting all these results, formulae of spectral indices are 
given in Table [T] and ID 


5. Applications 
5.1. The a ~ — 1 problem 


With the known expression of spectral indices, we would like to discuss the a ~ — 1 
problem hrst. From Table [H the fast cooling phase is the only possible case to explalin 
a ~ — 1. Since p > 2, a in the slow cooling “a” case is less than —3/2. Thus, we numerically 
investigate the pi value for different parameters. We define a new parameter, i.e.. 


^ U' 


2L. 


B 




( 20 ) 


which describes the ratio of the source photon energy density over the magnetic held energy 
density. We then plot pi as a function of r] for different values of ^ in Figure [TJ It is evident 
that the minimal value of pi decreases when ^ increases. This means that the Compton 
processes will strongly Fatten the electron spectrum in the low energy range. Note that the 
minimal value of pi can be 0 if £ goes to i nhnity . Correspondingly, the value of a in the band 
spectrum can reach —0.5. iDuran et al.l (120121) considered a in Equation flT5|) . and pointed 
out that the extreme value of a is —1. In our analysis, a ~ — 1, or equally pi ~ 1, is not a 
limit but a median value in the whole parameter spaces. The range of pi can be used to ht 
more GRB spectra. 


In order to show the parameter space more clearly, we plot the contour of pi in Figured 
where the x and y axes are logp and log.^, respectively. The green region shows a “boomerang” 
pattern in Figure [21 denoting for the pi ~ 1 zone. In the left wing of the boomerang, ^ can 
be several tens up to 10” {n > 4) when p is several. In the right wing, logp and log^ show 
a linear relation. We also want to know the value of Compton parameter of Yc in the same 
parameter space, since it denotes which radiation mechanism is dominant. The contour plot 
of Yc is given in Figure |3l One can observe that Yc is larger than 1 in the warm red and 
yellow color region, and is very small in the green and blue region. The separating zone 
between the warm and cold region is very close to the right wing of the boomerang, if Figure 
El and 121 are overlapped. So, pi can vary signihcantly when powers of the synchrotron and 
Compton scattering are comparable. 


One question arises naturally, if these parameter values rehect the true physical con- 
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Table 1. Spectral indices in Band function 


Spectral indices 

a 

/S 


Ep 

p 

<1 

Fast cooling 

pi+i 

2 

P2 + 1 

2 

P-1 

2 

hVm 

2^ + 1 

—2a — 3 

Slow cooling a 

p+1 

2 

P2 + 1 

2 

P2-P 

2 

huc 

-2a - 1 

-2A- 1 

Slow cooling b 

2 

3 

P+1 

2 

P 1 

2 6 

hUm 

2^+i 


Marginal case 

2 

3 

P2 + 1 

2 

P2 1 

2 6 

hl^rn^c 

2S-l-q 



Note. — The parameter S is defined as S = a — /3. In the “Slow cooling a” case, 
one has p < 3; while one has p > 3 in the “Slow cooling b” case. In the marginal 
case, one can only determine the sum of p and q. 


Table 2. Spectral indices in the double broken power law 


Spectral indices 

a 

/S 

7 

p q Ehi 

Eb2 

Fast cooling 

2 

3 

pl+1 

2 

P2 + 1 

2 

2{l3 - 7) + 1 -2f3 - 3 huc 

hUm 

Slow cooling 

2 

3 

p+1 

2 

P2 + 1 

2 

-2/3 - 1 2(/3 - 7) - 1 

hvc 


Note. — A double broken power law is possible in both the fast and slow cooling phases, p 
denotes the spectral index of injected electrons, q denotes the correction from the Compton 
processes. Ei,i and i?f ,2 are the broken energies in the double broken power law spectrum. 
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ditions in GRBs. We estimate the r] and ^ values for typical GRBs. rj has already been 
estimated to be 12 in the last section. Here, we choose to calculate rj = ■jpEp/MeC^ for 
reference, i.e.. 


V ~ 10(1 + ^)7p,3r2 


-1 


E„ 


0.5MeV 


( 21 ) 


The convention Q = 10”^Q „ is taken. The most common observed Ep is 300 keV flGoldstein et al. 
2 OI 2 I: iGruber et al.l 120141) . In the fast cooling phase, 7 ^ is the minimal Lorentz factor 7 ^^ of 
injected electrons, which is a free parameter unless some spe cihc accele ration mechanisms 
are specihed. For the shocked electrons, one has = OlO^eF (jSarilllOOSi) . Thus, rj is in the 
order of 10 for £e ~ 0.1. In the synchrotron radiation, the magnetic held B' can be estimated 
as 


B' ^ 0.4 X 10®(1 + 


Er, 


0.5MeV 


Gauss. 


( 22 ) 


Once the magnetic held B' is known, ^ is calculated from the dehnition, i.e.. 


«“l'lxl0-'(l + ^)-«*(s^) 


-2 


(23) 


It seems that £ is v ery small. This can be realized in the magnetic dominated jet model 
fiGhang et al. 2012b '). The way to increase ^ is quite limited. R is constrained by the variation 
time of GRBs. Ep and are observed variables. When we increase 7 ^, ^ is signihcantly 
enlarged. But, this also increases the value of 7 . One can avoid to increase rj by considering 
the low energy external photons, which is the case in the external Gompton model. In Figure 
m one observes that pi is smaller than 2 unless log^ is positive. Therefore, the left wing of 
the green zone is difficult to be located in the SSG model. 

In the external Gompton (EG) model, by increasing 7 ^ and decreasing in the same 

time, one can realize that rj is unchanged while ^ is greatly enhanced. When 7 ^ ~ 2 x 10® 
and ^source ~ 10 keV, one has r; ~ 10 and ^ ~ 100. The magnetic held is B' ~ 10 Gauss, 
which may be the shocked magnetic held in the jet. The resulted pi is in the green zone. 
Gorrespondingly, the Yc contour tells us that the power of Gompton is roughly the same 
with that of synchrotron, which only doubles the energy budget. By this constraint, we can 
discard the left arm of the boomerang, where the corresponding Yc is much larger than 1 . 
In some bursts, a black-body (B B) component plus the Band function can better ht the 
observed spectra IYu et al.l (120141 1. The BB bump in the spectra is found to be 10 keV. 
Gombining these information, the prompt emission of GRB can be explained by the reverse 
shock plus photosphere model. A reverse shock compresses the magnetic held and accelerates 
electrons in the jet, then electrons collide with the thermal photons and give synchrotron 
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radiation in the same time. This produces the observed Ep and low energy spectral index. 
The high value of 7 ^ rnay als o be due to particle acceleration in the magnetic reconnec¬ 
tion fjZhang &: Meszarosll2002l) . So, photospherical emission plus the magnetic reconn ection 
model is a also a hopeful theory to explain the a ~ — 1 problem fjZhang fc Yanll201ll) . 


The SSC model is not excluded completely, since the right wing of the boomerang is 
available. If 7 ^ increases to 10®, larger ^ and rj can be realized. In Equation ([23]), ^ depends 
on 'y'p. One observes that the pi = 1 contour is a straight line for logr^ > 2 . One needs a 
hne tuning of 7 ^ to land on the green zone. This strongly constrains all the relevant physical 
parameters. The relation between log.^ and logr; can also be described by a straight line. 
One hnds that log.^ = dlogr/ -|- b, where 6 ~ —11 denotes all other parameters except 7 ^ in 
Equations fET]) and fl25]l . We have marked this line with the white color in Figure |2J The 
crossing point of these two lines will determine the precise value of 7 ^. We hgure out roughly 
that Pi locates at 1 when 7 p 1.5 x 10^. The corresponding coordinate of this point is 
(log 7 ,log^) = (5.17,9.69) in Figure [T] The same point in Figure |3] corresponds to Yc ~ 1, 
which means that the power of the Compton process is the same with that of the synchrotron 
radiation. One also note that the pi ~ 1 region and Yd ~ 1 region agrees in a large area. 
Thus, we do not make an energy crisis to obtain pi = 1. 


Substituting 7 ^ 1.5 x 10^ back into Equations fl2l|) and fl2^ . one obtains that B' ^ 

10“^ Gauss and p ~ 1.5 x 10®. The Compton scattering processes are in the deep KN r egime. 
In the circum-stellar medium (CSM), B' is at the order of /i-Gauss I KnmarlI2009I) . The 
shocked magnetic held can be amplihed to hundred times larger. So, the derived magnetic 
held B' coincides with external shock model. With the price of 7 ), ~ 1.5 x 10^, the complex 
radiation mechanism can be a reasonable solution to the a ~ — 1 problem. However, the 
extreme high value of 7 ^ challenges the shock acceleration mechanism of electrons. The high 
energy electrons may be produced by the strong coupling between electrons and pro tons in 
the explosive jet, where the temperature is extremely high (jZhang fc Mfezarosll2002[) . Then 
these electrons are accelerated by the forward shock. A magnetic dominated outflow colliding 
with the CSM offers necessary gradients to explain all parameters obtained here, and needs 
future studies. 


We have used the typical parameters of GRB to calculate, but they are different from 
burst to burst, a is also different from burst to burst. The inclusion of Gompton processes 
offers a new method to fit the observed spectra. Recently, it was found that a ~ —0.7 for the 
time resolved spectra, which favors the synchrotron radiation model flYu et al.ll20141) . The 
light travel time effect, i.e. the observed photons are from electrons in different evolution 
stage, may also be a reason to explain a ~ — 1 , since —1 is about the averaged value of 
— 1/2 and —?)/2. Our analysis here also indicates that a can be any value between —2/3 and 
















14 


—3/2, if Compton processes are included. 


5.2. Phase transition 


In Table [H we have four possible phases available to explain the GRB spectra. Phase 
transitions between them may occur when 7 /^ and 7 ' vary with time. 7 /^ is determined by 
how electrons are injected, i.e. the acceleration processes. Except the slow cooling phase a, 
is related to the peak energy Ep. By the Equation fl20l) . the peak energy Ep of synchrotron 
origin is expresses as 


1 + Z MeC3/2 




(24) 


Here we include the redshift to count for the co smological effects. Ob viously, one observes 

1 /o 

that Ep oc , which is the Yonetoku relation (lYonetoku et al.l 1200411. If we change to 


the isotropic energy Eiso, this relation is the Amati rela tion flAmati et al.l 


Amati relation is evident by statistics of many bursts (lAmati et ah 


2002 


2008 


Yu et ah 


2008). The 


20141) . 


other parameters should be similar among different bursts. In the fast cooling and the slow 
cooling b phases, one can estimate by the following relation. 


dm 


= (1 + 


(25) 


In the “slow cooling a” case, Ep is related to 7 '. The critical Lorentz factor in Equation 
is calculated to be 

1 


7c = 




(26) 


which is obtained by requiring 
with energy Ec, i.e.. 


(Tt 1 + Yc 

The electrons with 7 ' radiate synchrotron photons 


E. = 9\/2- 


TT 


, hqeMeC^/‘^ 1 


(Tr 


1 + z^l + Yc 


)2^3/2^-3/2r6^_ 


(27) 


One observes that E^ oc E^ strongly depends on T, so it can vary fast during one 

burst. This relation is used to explain the anti-relation between the low energy peak and 


the b olometric luminosity in blazars, if the Ec is the observed low peak energy (iLvu et ah 

20141 ). 


Both y'^ and 7 ' can be changed during the prompt emission, y'^ signifies the injection of 
the photon energy, which marks the kinematic of the jet. It can be increased in the beginning 
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pulse of the burst. At the afterglow phase, no more energy is injected, 7 ^ will be decreased. 
Assuming 7 ' is invariant for a steady state, the variation of 7 ^ leads to the phase transition. 
Most probably, a slow to fast transition occurs in the flux arising phase, and a fast to slow 
phase transition happens in the flux decaying phase. 

The variation of 7 ' depends on the radiation power. Suppose that the system is in a 
slow cooling phase in the beginning, which means 7 c > 7 m- And the synchrotron radiation 
is in a steady state. As the photon density increases (the SSC case), the Compton scattering 
becomes important. The Compton scattering will increase the radiation efficiency, and 7 ' 
becomes smaller, see Equation fl29|) . When 7 ' is smaller than 7 ^,^, a slow to fast phase 
transition will occur. Electrons with 7 ^^ will quickly radiate away their energy, and cool 
down to the energy 7 '. Thus, the major population of electrons will accumulate at 7 '. Other 
cases of phase transition are also possible if and 7 ' both vary. All these phase transitions 
are indicated by the change of spectral indices and the peak energy. They can be used to 
explain the time-resolved spectra in many bursts. With the observed a, {3 and Ep, one can 
analyze the time resolved spectra of GRBs. 


As a concrete example, we consider the time resolved spactrum of GRB 080916G. GRB 
080916G is a lo ng burst well know n for its high redshift ~ 4.35 and extreme luminosity 
L-y ~ 10®'^ erg flAbdo et al.ll2009a|) . The minimal bulk L orentz factor was e stimated to be 
around 600 to meet the optical depth constraint < 1 flAbdo et al.ll2009al) . According to 
the variability time Af, the emission radius R is estimated to be around ~ 10^®cm. The 
peak energy evolves in time. In the hrst time interval a (0 ~ 3.58 s since the trigger time), 
Ep is about 440 keV. It goes up to 1.1 MeV in the second time interval (3.58 ~ 7.68 s), and 
then decreases. With these parameters, one can obtain that 


% « 1,2 X 10‘(1 + , 


7 ' ^ 3.4 X 10-3^ 


IMeV^ 

1 + ^ 7 , 54 ^ 2 -^ 16 - 


(28) 

(29) 


The minimal Lorentz factor 7 ^^ meets t he expect ed value of Lorentz factor for shocked 
electrons, i.e., blOSeE ~ 0.6 x 10 ^£e,-ir 2 fjSarilll998l ). The critical Lorentz factor 7 ' is much 
less than one. This is misleading since 7 c > 1 always holds. The derivation of 7 ' depends on 
the assumption ~ If t'dyrn 7c ^e much larger than one, or even larger 

than 7 p. The phase of the system is undetermined by just comparing the derived values. 

In the time interval a, a and (3 are observed to be —0.58 ± 0.04 and —2.63 ± 0.12, 
respectively. In the second tim e interval, a and (3 are —1.02 ± 0.02 and —2.21 ± 0.03, 
respectively flAbdo et al.l l2009al) . Then, the spectral indices do not vary signihcantly in 
the following time, and the peak energy decays. From Table [H if the burst is in the fast 
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cooling phase in the time interval a, one has p = 5.10 ± 0.32 and q = —1.84 ± 0.08. The 
spectral index of electrons is very large. Considering the slowing cooling b phase, one obtains 
that the spectral index of electrons is p = 4.26 ± 0.24. In the marginal phase, one obtains 
p + q = 3.26 ± 0.24. If Compton process is not signihcant, or equally g ~ 0, one has 
p = 3.26 ±0.24. In the second time interval, the emission zone are most probably in the fast 
cooling phase. With known a and /?, one obtains that p = 3.38 ±0.1 and q = —0.98 ± 0.04. 
Considering the continuously electron injection, common knowledge tells us that the spectral 
index of injected electrons should be roughly the same in different time intervals. The slow 
cooling case b is not preferred for this reason. So, the system is most probably in the marginal 
phase in time interval a. With constant p, the variation of q can account for the spectral 
indices in both the hrst and second time intervals. In both these two phases, Ep corresponds 
always to 7 ^. The revolution of Ep agrees with that of the time resolved luminosity. So, we 
have a high level conhdence to claim that a phase transition occurs in the prompt emission 
of GRB 080916C, i.e., transition from the marginal phase to the fast cooling phase. 


5.3. Time resolved spectra 


The time resolved spectra offers us more informat i on on the understanding of the ra¬ 
diation mechanism of the prompt emission. lYu et al.l (120141) presented the time resolved 
spectra of eight most energetic bursts with high signal-to-noise level. The distributions of 
spectral indices and Ep peak at —0.73, —2.13, and 374.4 keV, respectively. However, theses 
parameters are obtained from the total 299 spectra of eight bursts, the revolution of param¬ 
eters in each burst was not discussed. So, we aim to further analyse these parameters, and 
investigate whether the variation of spectral indices can be explained in the frame of complex 
radiation. We choose the data of three bursts for discussion (see Table A. 2 ., Table A.3., and 
Table A. 8 . in Yu et ah (2014) for reference), which are GRB 100724B, GRB 100826A, and 
GRB 130606B, respectively. The light curves of these three bursts are composed of sim¬ 
ple and clear pulses, and their Band spectra have better htting results. The htted spectra 
containing parameters with only upper limits are excluded. With these selection rules, the 
number of the selected spectra is 114 in total. These spectra are further classihed according 
to the pulse, we have seven pulses in total. All the pulses are illustrated in Table [3l We plot 
the /3 Ep^ CX Ep, and Ep ~ Ep relations for the three bursts, respectively. The peak flux 
Ep is calculated by Ep = EpN{Ep), see Equation flTSD . We use the linear relation y = a + bx 
to £t the these relations, and results are presented in notes of hgures. According to these 
relations, the phases of spectra will be analyzed in the following. 


GRB 100724B is composed of two pulses. Note that the second pulse is composed of 
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several relative small peaks in the light curve, but we still consider them to be in one pulse 
for convenience of classihcation. The plots in Figure H] indicate that /3 becomes smaller when 
Ep increases, while a is almost independent of Ep. The small slopes of two htting lines for 
(3 ^ Ep relation tell us that the Compton component may exist but not signify. The value of 
a also approaches —2/3 roughly. So, both pulses are of the marginal case. The coincidence 
of huc and hvm can be reflected by the Ep ~ Ep relations, which are not well linearly fitted. 
In the marginal phase, we can not determine p and q in principle, but we can obtain their 
constrains via the relation p + q = —2(3 — 2. The relation of p + g versus Ep is plotted in 
Figure [6l The linear httings are not good for the second pulse, this probably is due to our 
rough classification of pulses. 


GRB 100826A contains one big pulse and one small pulse. The plot oi (3 ^ Ep relation 
indicates that (3 decreases when Ep increases, and also a have the same trend. The Ep ~ Ep 
relation reveals that Ep is proportional to the EpN{Ep) flux, which is a sign of the fast 
cooling phase. Thus, both pulses of GRB 100826A are in the fast cooling phase. In the 
complex radiation mechanism, the difference between a and /?, i.e., S = a — (3 is invariant 
in principle. We plot the p ^ Ep relation in Figure IHl where p denotes the spectral index of 
injected electrons. It is evident that p remains almost invariant when Ep varies. The q Ep 
relation is also plotted in Figure |9l Except some exceptional points, there is a signihcant 
correlation between q and Ep. The intercepts of the htting lines are near —2. The two lines 
have different slopes, which indicates the intrinsic physical parameters are different for the 
two pulses. The correlation between q and Ep indicates that the Gompton scattering is in 
the deep KN regime. The different trend of p and q dependence on Ep can be considered as 
an evidence for the complex radiation mechanism. 


The light curves in GRB 130606B, see Fig.A.2. in (lYu et al.l 120141) . have four pulses in 
total. In Table [31 we combined the latter two pulses as one. This may bring some large errors 
in linear htting. In Figure flUl (3 and a have diherent dependence on Ep for three pulses. In 
pulse 1, the a and (3 are independent of Ep, and a approaches —2/3. The E Ep relation 
illustrates that pulse 1 is in the slow cooling b or marginal case. In the slow cooling b case, 
the injected electrons should have a spectral index p > 3, this phase can not be excluded. 
The nearly zero slope of the linear htting oi (3 ^ Ep relation indicates that the synchrotron 
radiation is dominant. In pulse b, both (3 and a show a slightly increasing trend when Ep 
increases. However, the linear htting has a large error. The most possible phase in pulse 2 
is the fast cooling with the pure synchrotron radiation, since the E Ep relation shows a 
nice linear htting, see Figure [TT] Pulse 3 have a large error correction for the (3 ^ Ep htting, 
while a has the similar behavior with pulse 1. The E Ep relation indicates a good linear 
htting with slope 1.636 ± 0.229. So, both slow cooling b and marginal case are possible in 
pulse 3. 
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The time resolved spectra show that prompt emissions of GRBs have variable phases 


from pulse to pulse. 


The Fp ~ 


Ep relations indicate that Ep mostly corresponds to hi>„ 


the slow cooling a phase is not evident in these three bursts. The spectra of GRB 100826A 
contains the feature of complex emission in a hig h signihcant manner. Such kind of feature is 
less evident in GRB 100724B and GRB 130606B. Frontera et ah ( 2012 1 found the correlation 
between a and Ep in GRB 980329, which is also a burst with single one pulse. One also 
observes that the flux rising phase and flux decaying phase have no differences in these 
relations during one pulse. This hints us that it is a better method to classify the time 
resolved spectra according to pulses. The Fp ~ Ep relations in their work als o show the nic e 


linea r httings in GRB 970111, GRB 980329, GRB 990123, and GRB 990510 flFrontera et ah 


201211 . Our results agree with theirs. The mar ginal phas e seems to be popular in the bursts. 


In this case, 7 ^ and 7 ' is close to each other. IYu et ahl (1201411 also revealed that there is a 
universal break ratio between and Uc, which is less than 10. Since 7 ^ corresponds to the 
acceleration processes, while 7 ' is related to the emission power of electrons, the agreement 
of them indicates the balance between the the acceleration and radiation processes. This is 
reflected by the fact that the Band function is suitable to describe the spectra in different 
time slices during one burst. 


6. Discussion and conclusion 

By considering both the SR and the Gompton processes, the continuity equation of 
electrons is investigated. We hnd that the spectrum of electrons is a broken power law in 
both the fast and slow cooling phases. The inclusion of the Gompton processes changes the 
spectral indices of photons. Analytical expression of spectral indices in four different phases 
of GRBs is given, which enriches the theoretical implications of the Band function. The 
a ~ — 1 problem can be solved in this new complex radiation frame. A detailed investigation 
of physical parameters shed lights on the prompt emissions of GRBs. In the EG case, the 
reverse shock plus the photosphere model is favored to explain the prompt emission. The 
nearly equal powers of Gompton and synchrotron radiations do not overburden the energy 
budget. Also, the SSG case can account for the spectral index and peak energy, if forward 
shock and strong coupling between electrons and heavy ions are considered. These arguments 
indicate that the popular shock models of GRBs can also answer the low energy spectral 
index problem. 

Phase transitions are also possible during the prompt emission of one GRB. We show 
that GRB 100826A experiences the phase transition, from the marginal phase to the fast 
cooling phase. We analyse the time resolved spectra of three bursts, which are GRB 100724B, 










- 19 - 


GRB 100826A and GRB 130606B, respectively. The plots of indices versus peak energy show 
many interesting features of the prompt emission. Different pulses have different correlation 
between indices and peak energies. Even in one burst, pulses can be in the different radiation 
phases. The correlations of /3 ~ Ep and a ^ Ep shows the signature of the complex radiation 
in the fast cooling phase, especially in GR B 100826A. The pure synchrotron radiation can not 
explain the revolution of spectral indices. iFrontera et al.l (l2012l) has reported the dependence 
of a on Ep in GRB 980329, which is also of one pulse burst. We guess that th e spectral indices 


depen dence on the peak energy can be found in the pulse dependent spectra. iTitarchuk et ah 


(120121 ) argued that the GRB spectra are formed two upscattering processes. The resulted 
spectral index of Gomptonization has relations with the bulk motion of the outflow and the 
otical d epth. Such model als o predicts a correlation between Ep and Liso for time resolved 
spectra flFrontera et al.l l2012h . However, (3 and a are not related in the Gomptonization 
theory. This can be considered as a special feature of the complex radiation in the fast 
cooling phase. 

In the complex radiation scenario, the derived p for the injected electrons is larger than 
3, e.g., p ~ 3.6 in GRB 100826A. This challenges the widely accepted electron acceleration 
mechanism. For shock acceleration, the typical va lu e of p is from 2.3 to 2.8 in the magnetized 
collisionless pair shock flSironi fc Spitkovskvll2009[) . iGuo et al.l (120141) have found that a hard 
power law of electrons can be formed in the relativistic magnetic reconnection, and p can 
be as hard as 1. The analytical way to obtain this hard spectra is also solved by the 
continuity equation in the unsteady state, where the injection, acceleration and escaping 
terms are considered. It was shown th at the hrst Ferm i acceleration can not produce a 
power law spectra without the injection (iGuo et ahl 120141) . These numerical results are not 
verified in observation, especially from GRB. The conflict here is a two sides thing. One 
side is that the soft injected spectrum questions the shock acceleration paradigm of GRB, 
specific physical configurations are needed to explain the steeper spectral index. Even in 
the pure synchrotron radiation, the derived p from Band fitting does not agree with the 
shock acceleration completely. The other side is that our analysis is questioned. In our 
derivation of the continuity equation, we consider only the single power law injection. If 
the injected electrons have a broken power law spectrum, the resulted photon spectrum is 
extended to double broken power law. Then, the Band function may be htted in the shock 
acceleration model. The more complicate spectrum can be obtained by complicating the 
injected spectrum. We did not consider the acceleration term and the time evolution term. 
If these terms are included, the power law solution still exists from analytical argument, and 
other corrections to the index can be expected. 


In our work, the Gompton induced spectra are not discussed. Roughly, the spectral 
index in the high energy range induced by Gompton processes are p 2 , since the electrons share 
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their energy with photons. The observation of the Compton indnced spectrnm is difhcnlt in 
GRBs, see Section 3. One can expect that this high energy spectral index depends on the 
synchrot ron peak energy and the IC peak energy. The Compton spectra are most evident in 
blazars. lAbdo et al.l (12 01 Obi ) fonnd a harder-brighter tendency in the ga mma ray band for 
certain specihc objects of the snbclasses of blazars, i.e., FSRQs and LBLs. ILvu et ahl (120141) 
indicated that FSRQs and LBLs are in the slow cooling phase, becanse there is an anti¬ 
correlation between Inminosity and peak energy. Combining them together, one can infer 
that the spectral index in the gamma ray band becomes softer when the peak energy increases 
for FSRQs and LBLs. For a samp le of bright blazars , a correlation between the spectral index 
and the peak energy was fonnd fjAbdo et ah 2010a ). The linear htting of the correlation is 


\ogVp^ = —4F -1-31.6 (F is the Fermi spectral index), see Figure 29 in fjAbdo et al.ll2010al) . 


The anti-correlation is against with that of GRB 100826A, which is obtained from the time 
resolved spectra. The disagreement does not make a real conflict, since blazars can be in 
the slow cooling phase, while GRBs are in the fast cooling ph ase. The investigatio n on the 


Gompton induced spectral indices will be given in the future f| Jiang et ak¬ 


in pre 


3)- 


One can conclude that the correlation between indices and peak energy (or luminosity) 
exists in both the time resolved and time integrated spectra, and this applies to both GRBs 
and blazars. Other correlations between spectral index and the peak energy are possible, see 
the time resolved spectra of the GRB 130606B. The spectral revolution is one main reason to 
cause the differences between the spectra of peak flux and the tim e integrated spec t ra, an d 
a large number of GRBs can have such behavior (see Figure 17 in iGoldstein et al.l f)2012l) ). 
Our analysis of the time resolved spectra in this work show the pulse dependent spectral 
revolution, and proposes a new manner to reveal the prompt emissions of GRBs. The 
richness of the complex radiation mechanism enables us to explain many observed spectral 
phenomena in astrophysical objects, especially the the GRBs and the blazars. 
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Pi values for different ^ 



In r| 


Fig. 1.— The index pi as a function of rj for ^ = 0.1,1,10,100 and 1000, respectively. 
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Fig. 2.— The contour plot of pi in the logr^ and log^ plane. 
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Fig. 3.— The contour plot of Compton Yq parameter. The numbers in the contour lines 
denote values of log Yc- 










24 


Table 3. The pulses of GRBs 


GRB 

Pulse 

Time interval (s) 

Number of spectra 

Phases 

100724B 

1 

-7.186 ~ 41.089 

10 

Marginal 


2 

41.089 ~ 130.458 

23 

Marginal 

100826A 

1 

-2.048 ~ 40.792 

36 

Fast cooling 


2 

40.792 ~ 98.549 

14 

Fast cooling 

130606B 

1 

-3.072 ~ 11.459 

7 

Slow cooling b/Marginal 


2 

11.459 ~ 26.218 

11 

Fast cooling 


3 

26.218 ~ 77.824 

13 

Slow cooling b/Marginal 


Note. — GRB pulses and their phases 


GRB 100724D GRB 100724B 




Fig. 4.— The plots of (3 and a versus Ep in GRB 100724B. In the left panel, the htting lines 
are described hy (3 = —1.491 ± 0.065 — (0.213 ± 0.086)i7p (black line) and f3 = —1.724 ± 
0.119 — (0.245 ± 0.313)i7p (red line) in pulse 1 and 2, respectively. In the right panel, the 
linear httings of a are given by a = —0.754 ± 0.062 — (0.078 ± 0.082)i7p (black line) and 
a = —0.710 ± 0.084 — (0.034 ± 0.222)Ep (red line), respectively. The value of Ep is given in 
unit of MeV in convention. 
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GRB 100724B 



Ep(MeV) 

Fig. 5.— The plot of Fp ~ Ep relation in GRB 100724B. The htting lines are given by 
logFp = -5.665 ± 0.197 + (1.078 ± 0.673)logEp (black) and logFp = -4.846 ± 0.181 + 
(2.583 ± 0.356)logi7p (red), respectively. 



26 


GRB 100724B 



Fig. 6.— The plot of p + g versus Ep relation in GRB 100724B. In pulse one, the linear 
htting is given hj p + q = 0.981 ± 0.130 + (0.426 ± 0.171)Ep, In pulse two, the htting line is 
given by p + q = 1.448 ± 0.237 + (0.489 ± 0.627)i7p. 
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GRD 100826A GRB 100826A 




Fig. 7.— The plot of jS and a versus Ep relation in GRB 100826A. In the left panel, the 
htting lines are described hj f3 = —1.692±0.074—(0.938±0.186)i7p (black) and [3 = —1.868± 
0.185 —(1.562±0.974)i7p (red), respectively. In the left panel, the linear fittings of a are given 
bya = -0.388±0.067-(0.670±0.168)Ep (black) and a = -0.379±0.116-(2.223±0.612)Ep 
(red), respectively. 
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GRB 100826A 



0.01 0.10 1.00 10.00 

Ep(MeV) 


Fig. 8.— The plot of Fp ~ Ep relation in GRB 100826A. The htting lines are given by logFp = 
-4.975 ± 0.077 + (1.686 ± 0.158)logEp and logFp = -5.471 ± 0.151 + (1.132 ± 0.162)logEp, 
respectively. 
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GRB 100826A GRB 100826A 




Fig. 9.— The plots of p and q versus Ep relation in GRB 100826A. In the left panel, the 
line httings of p versus Ep are described by p = 3.608 ± 0.196 + (0.536 ± 0.489)i?p and 
3.976 ±0.247— (1.327± 1.298)i7p in pulse 1 and 2, respectively. In the right panel, the linear 
htting of q are given hy q = —2.223 ± 0.135 ± (1.339 ± 0.335)i7p and q = —2.242 ± 0.233 ± 
(4.45 ± 1.224)i7p in pulse 1 and 2, respectively. 


GRB 130606B 


GRB 130606B 




Fig. 10.— The dependence of /3 and a on Ep in GRB 130606B. In the left panel, the 
three htting lines in pulse 1, 2, and 3 are /3 = —2.044 ± 0.044 — (0.056 ± 0.051)i7p, = 
-2.355 ± 0.360 ± (0.247 ± 0.517)Ep, and f3 = -2.55 ± 0.260 - (0.179 ± 1.21)Ep, respectively. 
In the right panel, the htting lines of a are given by —0.794 ± 0.068 — (0.223 ± 0.112)i7p, 
-1.32 ± 0.136 ± (0.257 ±0.196)Ep, and -0.572 ± 0.071 - (0.767 ±0.333)Ep in pulse 1, 2, and 
3, respectively. 
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GRB 130606B 
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Fig. 11.— The plot of flux on ~ Ep in GRB 130606B. The three htting lines are described by 
logFp = -4.990±0.144 + (1.216±0.401)logEp, logFp = -5.085±0.086 + (1.293±0.177)logEp, 
and logFp = —4.689 ± 0.183 + (1.636 ± 0.229)logi7p in pules 1, 2, and 3, respectively. 
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